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Spider capture silk is a biomaterial with both high strength and high elasticity, but the structural 
design principle underlying these remarkable properties is still unknown. It was revealed recently 
by atomic force microscopy that an exponential force-extension relationship holds both for capture 
silk mesostructures and for intact capture silk fibers [N. Becker et ai, Nat. Mater. 2, 278 (2003)]. 
In this Letter a simple hierarchical chain model was proposed to understand and reproduce this 
striking observation. In the hierarchical chain model, a polymer is composed of many structural 
motifs which organize into structural modules and supramodules in a hierarchical manner. Each 
module in this hierarchy has its own characteristic force. The repetitive patterns in the amino acid 
sequence of the major flagelliform protein of spider capture silk is in support of this model. 

[Phys. Rev. Lett. 94, 028104 (2005)] 
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The capture silk is a natural material produced by orb- 
web weaving spiders. Spiders rely on it to entrap flying 
preys H. Like the spider dragline silk and many other 
naturally occurring silks, the capture silk has a tensile 
strength that is comparable to steel; but, unlike steel, it 
is also extremely elastic, with the ability to be stretched 
to almost 10 times its relaxed contour length without 
breaking 0, 0] . This perfect combination of strength and 
extensibility conveys a high degree of toughness to the 
capture silk: its breakage energy per unit weight is more 
than 20 times that of a high-tensile steel 0- With the 
aim to produce synthetic silks with similar mechanical 
properties, materials scientists have devoted many exper- 
imental and computational efforts to the understanding 
of spider silk's structural organization 4J. Despite these 
painstaking efforts, the mechanism behind spider silk's 
remarkable strength and elasticity is still largely miss- 
ing, partly because of the difficulty to obtain high-quality 
crystallized structures of silk proteins. 

Single-molecule manipulation methods have been re- 
cently applied on spider silks to obtain very detailed in- 
formation on their force-extension response 13; • In a 
recent experiment, Hansma and co-workers '3| attached 
capture silk mesostructures (probably composed of a sin- 
gle protein molecule) or intact capture silk fibers to an 
atomic force microscopy tip and recorded the response of 
the samples to external stretching force. They found a re- 
markable exponential relationship between the extension 
X and the external force /, 

/cxexp(x/^), (1) 

where £, the length constant referred to in ^3], is a fitting 
parameter whose physical meaning needs to be decided 
(see below). £ — 110 ± 30 nm for a capture silk molecule 
and £ — 11 ± 3 mm for an intact capture silk fiber 0. 
The length constant of a silk fiber is about 10^ times 
that of a silk molecule; its relaxed contour length is also 



about 10^ times that of a molecule It appears that £ 
holds an approximately linear relationship with the con- 
tour length. 

The exponential force-extension curve is significantly 
different from the data observed duringstretching single 
double-stranded (ds) DNA molecules p] or single titin 
proteins 0, ■ The behavior of dsDNA can be under- 
stood by the wormlike chain model of entropic elasticity 
^3 J ^'^d that of titin by a two- level system coupled 
with entropic elasticity llj. Similar exponential force- 
extension data were also observed by Dessinges et al. 
when they stretched a single-stranded (ss) DNA molecule 
at low salt conditions (12i] . The data was explained as a 
result of the interpla y o f electrostatic repulsion and en- 
tropic elasticity [l2, Ujl- However, the success of this 
model depends on the specific ionic concentration (i.e., 
low salt and high electrostatic interaction); and it could 
not naturally reproduce the exponential stretching data 
at extremely high force (i.e., when the extension is larger 
than 1.1 times the ssDNA backbone length) if hi ghe r or- 
der deformation energy terms are not included [l^j- In 
the spider capture silk experiment, the exponential be- 
havior was observed at both fluid and air within a force 
range from about 10° piconewton (pN) to about 10^ pN 
3]. 

Equation IpJ indicates the following: (i) Because the 
capture silk is highly extensible, a great amount of ex- 
tra length must have been stored in its relaxed form, (ii) 
Since extension increases with force logarithmically, some 
fraction of the stored length must be easy to be pulled 
out, some fraction be harder to be pulled out, and till 
some other fraction be even harder to be pulled out. To 
model this kind of heuristic cascading responses, here we 
propose a hierarchical chain model for spider capture silk. 
In the hierarchical chain model, the polymer is composed 
of many basic structural motifs; these motifs are then or- 
ganized into a hierarchy, forming structural modules on 



more and more longer length scales. At the deepest hier- 
archy level hm, the structural motifs could be /3-s heets. 
/3-spirals, hehces or microcrystal structures 0. The 
interactions among some of these motifs are much more 
stronger than their interactions with other motifs, there- 
fore they form a structural module at the hierarchy level 
(/im— 1). These level- — modules are then merged 
into level-(/im— 2) modules through their mutual interac- 
tions. This merging process is continued; and finally at 
the global scale, the whole spider silk string is regarded 
as a single module of the hierarchy level h = 0. 

In nature, the structures of many biomaterials are in- 
deed hierarchically organized. As a composite material, 
chromosome is a mixture of DNA, histone, and other non- 
histone structural proteins 15]. The DNA molecule first 
wraps onto histone proteins to form nucleosome particles, 
the basic units of chromosome. With the help of HI his- 
tone, this linear sequence of nucleosomes are then coiled 
and folded to form the 30-nm chromatin fiber. With 
the help of other scaffold proteins, the chromatin fiber is 
then further coiled and folded at several levels to form the 
compact chromosome structure. As another example, the 
amino acid sequence of a protein first forms basic struc- 
tural modules of a- helix and /3-sheet, called secondary 
structures. By different ways of connections of these sec- 
ondary modules, the protein constructs various tertiary 
topologies that are critical for its specific biological func- 
tions. In a higher level, these tertiary domains then form 
quaternary structures of protein complexes consisting of 
multiple chains [l^ . Recent experimental and theoreti- 
cal studies have shown that the force-induced unfolding 
of protein molecules is indeed processed in a hierarchical 
way, i.e. the tertiary structure precedes the secondary 
structures to be pulled over at increasing external forces 
0, ^3 ■ Our hierarchical chain model may also serve as 
a framework to understand the mechanical property of 
these biomaterial systems. 

Consider a polymer of contour length Lq. It is regarded 
as the module of hierarchy level h — (the module Mq). 
Mq is composed of a tandem sequence of tuq subunits 
Ml of contour length Li = Lo/mo (Fig. P). Under the 
action of an external force field /, the extension of Mq 
is denoted as xo{f). It was observed that water-induced 
protein mobility is a significant contribution to capture 
silk elasticity 18]. Therefore, we decompose xq into two 
parts: (i) the extension Axq caused by the displacement 
and position-rearrangement of these mo subunits and, 
(m) the total extension niQXi caused by the inherent de- 
formations of these mo subunits: 



xoif) = Axo(/) + moXi{f). 



(2) 



Similarly, at the hierarchy level h — 1, each unit is 
itself composed of mi level-2 subunits M2 of length 
L2 = Lo/{mQmi). Therefore, xi{f) can be written as 
xi{f) — Axi{f) + miX2{f), where Aa;i(/) is the contri- 
bution to the extension of a level- 1 module due to the 




FIG. 1: The hierarchical chain model. At each hierarchy level 
h a structural module Mh is composed of a tandem sequence 
of submodules M^+i of hierarchy level h + 1. Under exter- 
nal stretching, Mh responds by (z) adjusting the arrangements 
of its rrih subunits and making them more aligned along the 
force direction, and (ii) extending these nih subunits. The 
thick broken lines between submodules of each hierarchy level 
indicate the existence of sacrificial bonds. 



displacement and position-rearrangement of its mi sub- 
units; and X2{f) denotes the extension because of the 
inherent deformation of a level-2 subunit M2. As this 
structural hierarchy is continued, we arrive at the follow- 
ing expression concerning the total extension x{f): 

momi . . .miiAxh+i{f)- (3) 

h=0 

It is reasonable that the average extension Axh{f) of 
a Mh contributed by the reorientation or rearrangement 
of its level-(ft,-|- 1) subunits is proportional to the relaxed 
contour length Lh of this module. To facilitate the fol- 
lowing analytical calculation, let us assume Axh{f) has 
the following non-linear form (we will show that the final 
force-extension relationship is not sensitive to the specific 
assumption made here): 



aLhf/fh, 



f <fh 
f>fh 



(4) 



where a is a dimensionless proportional constant; fh is 
the characteristic force needed to displace and rearrange 
the positions of the rrih submodules contained in a Mh 
(during this process some sacrificial bonds are broken). 
Equation is in agreement with the experimental ob- 
servation that, between adjacent rapture events, a 
capture silk responds to external stretching in a linear 
way. Denote AEh as the energy cost of breaking all the 
sacrificial bonds between a level-/i module's m^ subunits. 
From Eq. (g)) we know that fh+i/fh = nihAEh+i/ AEh- 
Consider a level-(h-|-2) module M^^j; if is in -^^;?+i which 
in turn is in M^. M^^^ feels an internal energy e due to 
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its interaction with other subunits in and it feels 

an external energy e' due to its interaction with other sub- 
units in but not in M^^^. Based on Fig.^ we know 
that AEfi^i = m/j+ie/2 and AE^ = mhrrih+ie' /2. The 
hierarchical organization of the polymer requires that 
e > e', so as to ensure that structural modules of shorter 
length scales will be formed earlier. Based on these con- 
siderations, we arrive at the following self-similar scaling 
form: 

h+i = {e/e')h = Ph (/3 = e/e' > 1). (5) 

The parameter (3 characterizes the degree of coherence 
in the modular organization of the polymer: a large /3 
value means that a submodule has much stronger internal 
interactions compared with its external interactions. 
From Eqs. and Q we find that when fh-i < f < fh 

Equation © therefore recovers the experimental expo- 
nential force-extension relationship of Eq. |^ with 



The length constant i is proportional to the relaxed con- 
tour length Lq of the whole polymer, in consistence with 
Ref. 0. 

Figure 12 demonstrates the numerically calculated 
force-extension curve based on Eqs. iQ, and ©. As 
a comparison, the experimental data (3| on intact spider 
capture silk is also shown. As /? ~ 2 and the exprimen- 
tal exponential force range is roughly from 6 x lO^^A^ 
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FIG. 2: Exponential force-extension relationship for the hi- 
erarchical chain model. Equation ||1J is used in the numerical 
calculation. The parameters are /o = 10~* N, a = 0.3, and 
P = 2 (the upper curve) or /3 = 1.75 (the lower curve). Ex- 
tension is in units of Lq. Symbols are experimental data from 
Fig. 4 of 0. 



to 10~^iV, it appears that 4-5 levels of hierarchy were 
probed. 

The exponential relationship shown in the figure is in- 
sensitive to our particular assumption Eq. as long 
as the elastic response at each hierarchy level is non- 
linear and bounded. As an example, the solid curves 
in Fig. 121 show the resulting force-extension relationship 
when Eq. is replaced by 

Axhif) = - exp(-//A)). (8) 

The same exponential behavior as in Fig. |2is obtained. 
However, the hierarchical scaling form Eq. |SJ| is needed 
for the exponential force-extension correlations. For ex- 
ample, when Eq. (|SJ) is replaced by a power-law, fh oc 
fah'^ , the response is not exponential (the dotted line in 
Fig.O. We also noticed that, when in Eq. |SJ the param- 
eter (3 is not a constant but fluctuates over some finite 
range of /? > 1, the resulting force-extension curve is still 
exponential (Fig.|3| dashed lines). 

In summary, we have developed a hierarchical chain 
model to understand the strength and elasticity of spider 
silks. Remarkably, this simple model was able to repro- 
duce the peculiar exponential force-extension response of 
spider capture silk reported by Becker et al. The 
model can also be used as a framework to understand 
the elasticity of other spider silks and other biopolymers 
with hierarchically organized structures. 
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FIG. 3: The force-extension relationship of a hierarchical 
chain is insensitive to the assumption made to the response 
Axh{,f) of Eq. Q, but is sensitive to the hierarchical scal- 
ing form of the characteristic force fh- The solid lines are 
obtained by assuming Axh{f) has the form of Eq. (|SJ, while 
other parameters are the same as those in Fig. |5| The dot- 
ted line shows the change in the force-extension curve when 
additionally a power-law form of fh = foh'' with 7 = 3.0 is 
assumed for fh- The dashed lines are obtained by assum- 
ing Eq. (|HJ and Eq. (|KJ, with (3 fluctuating uniformly within 
[1.25, 2.25] (the lower curve) and within [1.5, 2.5] (the upper 
curve) . 
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Becker et al. |3] have proposed an alternate and inter- 
esting idea to model the spider silk as a branched network 
of interconnected springs. In their model, the system re- 
sponds to external stretching by first unfolding the sin- 
gle spring at the root level, then the m springs at the 
first branching level, then the m? springs at the second 
branching level, and so on. The hierarchical chain model 
developed here is different from Becker et al.'s model. 
First, the molecule in the hierarchical chain model con- 
sists of a tandem sequence of structural modules of differ- 
ent length scales. Therefore, no assumption of branching 
structure is made in our model. Second, the response 
of the chain to external perturbations is in a hierarchi- 
cal manner. If the external force is small, only those 
structural units of length scale comparable to the whole 
polymer length will be displaced and rearranged; struc- 
tural units at short and moderate length scales will re- 
main intact. As the external perturbation is elevated, 
additional structural units at more shorter length scales 
are also deformed. Through such a hierarchical orga- 
nization, a single polymer chain can respond to a great 
variety of external conditions; at the same time, it is able 
to keep its degree of structural integrity as high as possi- 
ble. This hierarchical modular structure also indicates a 
broad spectrum of relaxation times. The modules at the 
shorter length scales will have much shorter relaxation 
times and will be refolded first when the external force 
decreases. This gap in relaxation times ensures that, af- 
ter extension, the spider capture silk will return to its 
relaxed state gradually and slowly. This is a desirable 
feature for spider capture silk, because a too rapid con- 
tract following the insect's impact would propel the vic- 
tim away from the web. 

The simple hierarchical chain model, while appealing, 
needs further experimental validation. This model seems 
to be supported by recent genetic sequencing efforts. By 
analyzing the cDNA sequence of the major protein of spi- 
der capture silk, the flagelliform protein, it was revealed 
that the amino-acid sequence of flagelliform has a hier- 
archy of modularity jjlj . At the basic level, the 
flagelliform sequence is consisted of three repetitive mod- 
ules (motifs): (i) the GPGGX motif (Gly-Pro-Gly-Gly- 
X, X e {Ala,Ser,Tyr,Val}); {ii) the GGX (Gly-Gly-X, X 
G {Ala,Ser,The}); {Hi) the highly conserved spacer motif 
of length 28 amino-acids. At the next level, an ensem- 
ble motif is formed, which is a tandem array of 43—63 
GPGGX followed by 6—12 GGX, the spacer and another 
GGX At the even higher level, the ensemble motif 
then repeats itself about 14 times to form the flagelliform 
monomer. The variable residues X are not randomly dis- 
tributed along the protein sequence (l5] . which may en- 
code important structural information. The structures 
of spider capture silks (and other spider silks) therefore 
have the potential to be hierarchically organized. 

For spider capture silk, one important experiment will 
be to check the validity of Eq. ^ by single-molecule 



force spectroscopy. The characteristic forces ft niay cor- 
respond to the forces at the saw-tooth rapture events ob- 
served by Hansma and co-workers 0. Although the ex- 
perimental curves of Ref. is in agreement with Eq. 
more systematic investigations are necessary to draw a 
solid conclusion. 



We are grateful to H. Hansma for communicating 
her experimental results with us prior to publication 
and for the experimental data used in Fig. [51 and to 
D. E. Makarov for comments on the manuscript. Thanks 
are also due to R. Lipowsky, Z. C. Ou-Yang, and J. Skol- 
nick for their support. 



[1] F. VoUrath, Sci. Am. 266, 70 (1992). 

[2] J. M. Gosline, P. A. Guerette, C. S. Ortlepp, and K. N. 

Savage, J. Exp. Biol. 202, 3295 (1999). 
[3] N. Becker, E. Oroudjev, S. Mutz, J. P. Cleveland, P. K. 

Hansma, C. Y. Hayashi, D. E. Makarov, and H. G. 

Hansma, Nature Materials 2, 278 (2003). 
[4] D. Kaplan, ed., Silk Polymers: Materials Science and 

Biotecfinology, ACS Symposium Series (Am. Chem. Soc, 

Washington DC, 1994). 
[5] E. Oroudjev, J. Soares, S. Arcidiacono, J. B. Thompson, 

S. A. Fossey, and H. G. Hansma, Proc. Natl. Acad. Sci. 

USA 99, 6460 (2002). 
[6] S. B. Smith, L. Finzi, and C. Bustamante, Science 258, 

1122 (1992). 

[7] M. Rief, M. Gautel, F. Oesterhelt, J. M. Fernandez, and 

H. E. Gaub, Science 276, 1109 (1997). 
[8] K. S. Z. Kellermayer, S. B. Smith, H. L. Granzier, and 

C. Bustamante, Science 276, 1112 (1997). 
[9] J. F. Marko and E. D. Siggia, Macromolecules 28, 8759 

(1995). 

[10] H. Zhou, Y. Zhang, and Z.-C. Ou-Yang, Phys. Rev. Lett. 

82, 4560 (1999). 
[11] M. Rief, J. M. Fernandez, and H. E. Gaub, Phys. Rev. 

Lett. 81, 4764 (1998). 
[12] M. N. Dessinges, B. Maier, Y. Zhang, M. Peliti, D. Ben- 

simon, and V. Croquette, Phys. Rev. Lett. 89, 248102 

(2002). 

[13] Y. Zhang, H. Zhou, and Z.-C. Ou-Yang, Biophys. J. 81, 
1133 (2001). 

[14] F. Vollrath and D. P. Knight, Nature 410, 541 (2001). 
[15] K. E. van Holde, Cfiromatin ( Springer- Verlag, New York, 
1989). 

[16] C. Branden and J. Tooze, Introduction to Protein Struc- 
ture (Garland Pub., New York, 1999), 2nd ed. 

[17] H. Lu, B. Isralewitz, A. Krammer, V. Vogel, and 
K. Schulten, Biophys. J. 75, 662 (1998). 

[18] K. M. Bonthrone, F. VoUrath, B. K. Hunter, and J. K. M. 
Sanders, Proc. R. Soc. Lond. B 248, 141 (1992). 

[19] C. Y. Hayashi and R. V. Lewis, J. Mol. Biol. 275, 773 
(1998). 

[20] P. A. Guerette, D. G. Ginzinger, B. H. F. Weber, and 

J. M. Gosline, Science 272, 112 (1996). 
[21] C. Y. Hayashi and R. V. Lewis, Science 287, 1477 (2000). 



